Nonlinear viscosity and velocity distribution function in a simple longitudinal flow 
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A compressible flow characterized by a velocity field u x (x, i) = ax /(I + at) is analyzed by means 
of the Boltzmann equation and the Bhatnagar-Gross-Krook kinetic model. The sign of the control 
parameter (the longitudinal deformation rate a) distinguishes between an expansion (a > 0) and 
a condensation (a < 0) phenomenon. The temperature is a decreasing function of time in the 
■ former case, while it is an increasing function in the latter. The non-Newtonian behavior of the 

gas is described by a dimensionless nonlinear viscosity n* (a*), that depends on the dimensionless 
longitudinal rate a*. The Chapman- Enskog expansion of rf in powers of a* is seen to be only 
asymptotic (except in the case of Maxwell molecules). The velocity distribution function is also 
studied. At any value of a* , it exhibits an algebraic high- velocity tail that is responsible for the 
divergence of velocity moments. For sufficiently negative a*, moments of degree four and higher 
may diverge, while for positive a* the divergence occurs in moments of degree equal to or larger 
than eight. 

PACS numbers: 05.20.Dd, 47.50. +d, 05.60.-k, 51.10.+y 



I. INTRODUCTION 
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One of the most challenging problems in nonequilibrium statistical mechanics is the understanding of transport 
' ^ ■ properties in fluids beyond the scope of the Navier-Stokes (NS) constitutive equations. As part of the NS constitutive 
equations, Newton's law establishes a linear relationship between the irreversible momentum flux and the velocity 
gradients, namely 



+ d V ' u ^j -CnsV-u% (i) 

where P^- is the pressure tensor, p = (l/d)TrP is the hydrostatic pressure, d > 2 being the dimensionality of the 
| system, u is the flow velocity, ?/ns is the shear viscosity, and Cns is the bulk viscosity jj| . The linear law (0) only holds 
£SJ . for small hydrodynamic gradients, i.e., when the typical distances over which the hydrodynamic quantities change are 
much larger than a characteristic microscopic length (such as the mean free path in the case of gases). Otherwise, 
Eq. (H) no longer holds, a situation usually characterized by the introduction of a generalized or nonlinear viscosity 
that depends on the hydrodynamic gradients H . 

The nonlinear viscosity has been extensively studied in the so-called uniform shear flow, which is characterized by 
a linear velocity field u x = ay and uniform density and temperature (^-Q). This is an example of an incompressible 
flow B, since V • u = 0. Recently, some attention has been devoted to viscous longitudinal flows of the form 
u(r, t) = u x (x,t)x |]§|-|l2||. The simplest example of such compressible flows is characterized by a linear velocity 
profile, i.e., u x (x, t) = A(t)x, and uniform density n and pressure tensor P |T^,[l4|]. In that case, the balance equations 
Q , for mass and momentum read 



dn . d(nA) „ ., 

whose solution is 

W-lhi- = T+ai' < 3 > 

where a is an arbitrary constant that represents the (initial) longitudinal deformation rate and no > is the initial 
density. In this case, Newton's law ([!]) becomes 

+ J A. (4) 
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This simple flow is known as homo- energetic extension and, along with the uniform shear flow, is a particular case of 
a more general class of homo-energetic affine flows characterized by d 2 Ui/ dxjdxu = In the flow defined by Eqs. 
(||) the longitudinal deformation rate a is the only control parameter determining the departure of the fluid from its 
equilibrium state, thus playing a role similar to that of the shear rate in the uniform shear flow state. On the other 
hand, in contrast to the uniform shear flow, the sign of a plays a relevant role and defines two distinct situations. The 
case a > corresponds to a progressively more slowly expansion of the gas from the plane x = into all of space. 
Given a layer of width 5, the flux of particles leaving the layer exceeds the flux of incoming particles by naS/(l + at) 
and, as a consequence, the number of particles inside the layer decreases monotonically with time. As time progresses, 
the system becomes more and more rarefied until no particles are left in the long-time limit, i.e., lim^oo n(t) = 0. 
On the other hand, the case a < corresponds to a progressively more rapidly condensation of the gas towards the 
plane x = 0. The latter takes place over a finite time period t = \a\ . However, since the collision frequency rapidly 
increases with time, the finite period t = | a | 1 comprises an infinite number of collisions per particle (see below). 

Equations (||) apply regardless of the initial density of the fluid. On the other hand, a kinetic description based on 
the Boltzmann equation is valid only for densities such that na d is much smaller than I, where a is a characteristic 
distance measuring the effective size of the molecules. Let us call n a ~ a~ d a characteristic density beyond which 
noticeable deviations from the Boltzmann equation can be expected. Thus, even if no <C n a , there exists a finite time 
t a = | a | 1 ( 1 — no/n a ) beyond which the Boltzmann description ceases to be applicable in the case a < 0. This time 
t(j can be made arbitrarily close to the maximum time |a| _1 by formally taking the limit no/n a — > 0. 

The aim of this paper is to carry out a detailed and self-contained analysis of the noncquilibrium behavior of a 
dilute gas under the longitudinal flow characterized by Eqs. (||), for arbitrary sign and magnitude of the control 
parameter a. The study is performed by using the tools of kinetic theory, namely the Boltzmann equation and the 
Bhatnagar-Gross-Krook (BGK) kinetic model, and deals with the nonlinear viscosity, as well as with more general 
velocity moments and the velocity distribution function. Most of the results are derived for arbitrary dimensionality 
and for the general class of repulsive potentials of the form r"^ with u > 2(d — 1). The Boltzmann equation for the 
problem is considered in Sec. O. Since the density is known, cf. Eqs. (g), one can focus on the probability distribution 
of velocities. In addition, the distribution becomes uniform when the velocities are referred to the local Lagrangian 
frame moving with the flow velocity u x (x, t). As a consequence, the original problem can be mapped onto that of a 
uniform system with a stationary density and subject to the action of a non-conservative driving force; also, there is a 
nonlinear relationship between the time variables in the original and the equivalent systems. To proceed further, the 
Maxwell interaction, /i = 2(d — 1), is assumed and the time evolution of the pressure tensor is exactly obtained. The 
long-time behavior allows one to identify the nonlinear viscosity as a function of the longitudinal deformation rate. 
When the velocities are scaled with the (time-dependent) thermal velocity, the distribution obeys in the long-time 
limit a steady-state Boltzmann equation with the addition of a second non-conservative force playing the role of a 
thermostat. The exact fourth-degree velocity moments are then derived as functions of the longitudinal rate and are 
seen to diverge in the case of condensation for states sufficiently far from equilibrium. The picture is complemented in 
Sec. Ill by the solution of the BGK kinetic model. In the case of Maxwell molecules, the distribution function exhibits 
an algebraic high-velocity tail that is responsible for the divergence of the moments. The solution predicts that 
moments of sufficiently high degree can also diverge in the case of expansion. In addition, the distribution function 
becomes infinite in the limit of zero velocity in far from equilibrium states. Section III also analyzes the nonlinear 
viscosity for more general repulsive potentials, including hard spheres. Evidence is given about the non-convergent 
(but asymptotic) character of the Chapman-Enskog expansion, Maxwell molecules being an exception. On the other 
hand, the nonlinear dependence of the generalized viscosity is practically insensitive to the interaction potential in the 
case of positive deformation rates. For negative rates, however, the influence of the potential is not small, especially 



near the maximum value of the viscosity. The paper ends with a summary of the main conclusions in Sec. IV 



II. BOLTZMANN DESCRIPTION FOR MAXWELL MOLECULES 



Let us consider a dilute gas subject to the homo-energetic extension flow described in the previous Section. All 
the relevant information is contained in the one-particle velocity distribution function f(x,v,t). In particular, the 
number density n, the flow velocity u, and the pressure tensor P are obtained in terms of velocity moments of /, 

n= [ dv/, rai= / dvv/, P = m / dvVV/. (5) 



Here m is the mass of a particle and V = v — u is the peculiar velocity. The trace of P gives the hydrostatic pressure 
p = (I/<i)TrP, which is related to the temperature T through the equation of state p = nksT, fee being the Boltzmann 
constant. The time evolution of / is governed by the Boltzmann equation Jig] 
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i+v*l)f = JlfJl («) 

where J[f, /] is the nonlinear Boltzmann collision operator, whose explicit expression will be omitted here. 

As happens in the uniform shear flow p|,|l6[|, the velocity distribution function f(x,v,t) becomes spatially uniform 
when the velocities are referred to a Lagrangian frame moving with the flow, i.e., f(x,v,t) — /(V, t), where V = 
v — u(x, t) is the peculiar velocity. The Boltzmann equation (0) for this flow can then be written as 



8t 8V X 



where 



/(V,T) = -^/(z,v,i), r = a- 1 ln(l + at). (8) 
n(t) 

Equation (]?]) can be interpreted as corresponding to a homogeneous gas with a velocity distribution / and subject 
to the action of a non-conservative force —maV x x. Note that the density and pressure tensor associated with / are 
n = no and Py = (no/n)Pij, respectively. In fact, / is proportional to the probability distribution of velocities. 
The time variable r = f Q dt'n(t')/no is a nonlinear measure of time scaled with the number density; this variable is 
unbounded even if a < 0, since in that case r — > oo when t — > | cz | 1 . It must be emphasized that Eqs. (p|) and (^) are 
fully equivalent in the present problem. On the other hand, Eq. ^ has the advantage of describing a uniform system 
with a constant density. The prize to be paid is the introduction of a driving force that acts as a longitudinal "drag" 
force in the case of expansion (a > 0) and as a "pushing" force in the case of condensation (a < 0). Of course, every 
solution to Eq. (Q) can be mapped onto a corresponding solution to Eq. (||). 

The Boltzmann equation (Q) cannot be solved by analytical tools in general. On the other hand, its associated 
hierarchy of moment equations can be recursively solved in the special case of Maxwell molecules, i.e., particles 
interacting via a repulsive potential of the form <p(r) oc r - 2 ( d - 1 ) j i n which case the collision rate is independent of the 
relative velocity of the colliding particles [ ]l7| j. In particular, Eq. (Q) yields a closed set of equations for the elements 
of the pressure tensor in the case of Maxwell molecules, namely 

d _ 2a ~ 

^P + -rPxx = 0, 9 
or a 



dr x 



-v (P xx -p), (10) 



where p = n^ksT = (l/d)TrP and vo oc no is a constant that plays the role of an effective (initial) collision frequency. 
More explicitly, Vq — p/^ns- Equation (||) is not but a condition expressing the conservation of energy. As for Eq. 
(|l0|), it must be emphasized that it is exact for Maxwell molecules in our problem, i.e., no approximate truncation 
scheme (such as Grad's moment method |0|) has been applied. This is a consequence of the fact that the collisional 
velocity moments of a certain degree do not involve moments of a higher degree in the case of the Maxwell interaction 



JL3|. The right-hand-side of Eq. (10) represents the (bilinear) collisional moment of 

In this case of Maxwell molecules the time variable r is just proportional to the average number of collisions per 
particle between and t, namely r = v^ 1 f* dt'v(t'), where v = p/^ns = {n/no)VQ is the time-dependent collision 
frequency. Thus, as said in Sec. Q, every particle experiences an infinite total number of collisions between the initial 
time and the finite interval t = | a | 1 when a < 0. After many collision times (r 3> Vq 1 ) both p and P xx behave as 
p, P xx ~ exp[— \(cl*)vot\, where a* = a/fo = A{t)/v(t) is the reduced longitudinal rate and A(a*) is the smallest root 
of the quadratic equation 

A 2 - (2a* + 1)A + ^a* = 0, (11) 



A(a*)=a* + i-v(" +-) --,<>■ ('->) 
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The second root is obtained from (|lj) by changing the sign of the radical and is only relevant in the transient stage. 
Consequently, 

ft, = dXia*) 

t^oo p 2 a* v ; 

In terms of the real time, one has an algebraic behavior for the temperature, T(t) ~ (1 + at)~ x ( a '/ a . It is important 
to note that the sign of A is the same as that of a*. This means that the temperature monotonically decreases in 
time if a* > (expansion, drag force), while it increases if a* < (condensation, pushing force). In addition, A 
is a monotonically increasing function of a* that behaves as A ~ 2a* + (d — l)/d in the limit a* — > — oo and as 
A w [1 - (d - l)/2da*]/d in the limit a* -> oo. 

Based upon Eq. A), we define in the long-time limit a (dimensionless) nonlinear viscosity as 



P - ftz 



2(d-l) ^ NS 



1 - — I , (14) 



2(d- l)a* \ p 

where we have taken into account that the bulk viscosity Cns vanishes in a low-density gas |l5| . Using Eq. ([l3"|), we 
simply get 

d 2a* - d\(a*) 

" ( ° )= d-l 4a* 2 • (15) 

In the three-dimensional case (d = 3), this result coincides with the one derived by Karlin et al. [|1 for Maxwell 
molecules by applying the invariance principle under the microscopic and macroscopic dynamics in the context of 
Grad's method. It should also be noted that, although by a somewhat different route, most of the above results (for 
d = 3) were first derived by Galkin more than thirty years ago JT^Jl4]| . 

The behavior of the nonlinear viscosity ( |l5| ) for small longitudinal rates is if ss 1 — [2(d — 2)/d]a*. More generally, 
?7* can be expressed as a series expansion in powers of a*: 

oo 

V *(a*) = l + Y / Cna* n . (16) 



This is just the specialization of the Chapman-Enskog expansion |15| , |18| to the simple viscous longitudinal flow. In the 
case of Maxwell molecules, a* = is a regular point of 77*, so the expansion ( p^) is c onvergent, although with a finite 
radius of convergence (\a*\ < |) due to a branch point at a* = ~(d—2)/2d±iy / d — 1/d. The knowledge of the explicit 
expression of 77*, Eq. (fl5l), allows one to get also its asymptotic behaviors for large |a*|; they are n* ss [d/2(d— l)]a* _1 
for a* > and 77* ss (d/2)\a*\~ 1 for a* < 0. This implies that P xx /p — * when a* — * +00, while P xx /p — > d when 
a* — > —00. In the former limit all the particles tend to move along the transverse directions, while in the latter they 
tend to move along the longitudinal direction. The shape of n*(a*) for d = 3 in the range — 2 < a* < 2 is shown in 
the next Section (cf. Fig. ||). 

Because of the symmetry of the problem, one expects that the heat flux is an irrelevant quantity that, even if it 
is initially different from zero, asymptotically decays in the long-time limit. Let us analyze this point more carefully. 
Taking third-degree moments in Eq. (^), we get 

^-q x + a(q x + 2M xxx "j = -~2/ §x, (17) 



d ( 3 1 

-g^M xxx + 3aM xxx = —vq I —M xxx — —q x ) , (18) 



where 



m I ._ n ~ — m 



j J dv r-'r,/. M xra .iv 1;/. (19) 

In Eqs. ( ]l7| ) and ( ]l8| ) use has been made of the third-degree collisional moments for the three-dimensional Maxwell 
interaction [ p^[ . For long collision times (r 3> ^o^ 1 ), the third-degree moments behave as q x ,M xxx ~ exp[— o;(a*)i/oT], 
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where oj(a*) = 2a* + y§ — y (a* + ^) - a*. If a* < — § + ^ —0.271, to is negative and then the heat flux grows 
in time. Apparently, this seems to contradict our expectation about the irrelevance of the heat flux in our problem. 
The solution to this paradox lies in the fact that also the temperature grows if a* < 0; indeed, what is relevant 
is not the absolute value of the heat flux but its value relative to the third power of the thermal velocity, namely 
q x /mn (2kBT/m) 3 / 2 ~ exp{— [uj(a*) — §A(a*)]^ T}. The difference uj(a*) — §A(a*) is always positive, but goes to 
zero as j^|a*| _1 in the limit a* — > — oo. In general, the smaller the value of a* the longer the transient period before 
the heat flux, conveniently scaled with the thermal velocity, has decayed to zero. 

In the long-time limit, not only the reduced elements of the pressure tensor Pij/p reach well-defined stationary values 
that depend on the reduced longitudinal rate a* , but the same happens with the distribution function when properly 
nondimensionalized with the temperature. This is just a statement on the validity of the "normal" or hydrodynamic 
regime, that applies for sufficiently long times []l5 18 1. To be more precise, let us introduce a reduced velocity £ and 
a reduced distribution function $ as 



1/2 



V, $(£; a*) = n 1 lim 



2k B T(r) 



(1/2 



/(V,r). 



2k B T(r)_ 

In this hydrodynamic regime, the Boltzmann equation (0) becomes (for Maxwell molecules) 



dt; x 



X(a*) d 
2 di 



(20) 



(21) 



As before, the first term on the left side represents a driving force. The second term can be interpreted as a thermostat 
force (HQ] that compensates for the heating (a* < 0) or cooling (a* > 0) effect produced by the former. Equation 
( pl| ) yields directly the "stationary" values of the (reduced) second-degree moments, namely 



1 



i 



2 1 + 2a* 



A 



(&.) = 



d-l 1 

2 1 - A' 



(22) 



where = £ 2 — £ 2 . The consistency condition (£ 2 ) = d/2 leads again to Eq. (|iT|). As noted before, (£ 2 ) — > in 
the limit a* — > +oo and — > in the opposite limit a* — > — oo. This means that all the particles move (in the 
Lagrangian frame) along directions perpendicular to the flow when a* — > +oo, while they move parallel to the flow 
direction when a* — > — oo. 

Now we can go further and consider the fourth-degree moments (£, i), (£x£±.)> an d (Cl)- Making use of the fourth- 
degree collisional moments for three-dimensional Maxwell molecules |13] , Eq. (pl|) gives rise to 



35 

54m-91 /<r2\2 



35 
3(36w- 



70 



35 

9(7-3iu) 
70 



(23) 



2(a*-A)<g£> = - 



144m+49 
210 



343- 162m /c2\2 
105 VSa 



+ 



r + 



Sw-7 I cA 
35 v*x 

469-216ui f j2 
140 v»x 



210 \ S -L/ 
3(36m-49) 



(24) 



o\/t4\_ 3«-+28 /,r4\ i 8(7-3m) /c 4\ , 4(18m-7) /c 2 c 2 \ 

- ^ A \?±; - — h log - \w H — log — 

i 2(154^-81) /C 2\2 18(7-3t») /C 2\ , 3(56-9t») / oc -n 
H 105 \W 35 ^*/ H 70 ' W 

where w ~ 1.8731 is the ratio of two eigenvalues of the linearized collision operator p9[ . The solution of this linear 
set of algebraic equations is 

(£) = 3 ^^y [-1296(7 - 3w)A 5 + 54(217 - 48w)A 4 - 9(327w + 602)A 3 

+33(129™ + 14)A 2 - 25(75™ - 14)A + 245w] , (26) 

1 — 3A 

= [ 43 2(7 ~ 3 W )A 4 - 54(64 w - 91)A 3 

+9(363w - 392)A 2 - 10(114™ - 35)A + 245w] , (27) 
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= — ^) [324(7 - 3w)X 4 - 90(6w + 7)A 3 

+9(143w + 28)A 2 - 10(114u; - 35)A + 245u>] , (28) 

where 

A(o*) = (1 - A) 2 [-432(7- 3w)A 4 - 54(8w - 7)A 3 

+9(159w - 56)A 2 - 10(114u; - 35)A + 245w] . (29) 

Upon writing Eqs. (^6|)-([29]) we have made use of ( p2| ) and have eliminated a* in favor of A. In the limit a* — > +oo, 
i.e., A — > |, we have (£ 4 ) = 9(56 — 9u;)/2(14io + 9) ~ 5.71 and (£ 4 ) = (? 2 £ 2 ) = 0, in agreement with a vanishing 
population of particles moving along the longitudinal direction. A more interesting situation occurs in the domain of 
negative a* (condensation case). As a* becomes more and more negative, the three moments (£ 2 £i), an d 
grow monotonically and eventually diverge when a* approaches a critical value a* ~ —1.599 (which corresponds to 
the root A c ~ —2.607 of the function A). This singular behavior of the fourth-degree moments also occurs in the case 
of uniform shear flow |20| , pl|| and is an indication of an algebraic high- velocity tail in the distribution function p2| , p3| . 
We will return to this point later on. The moments (£ 4 ), (QQ), and are plotted in Fig. [1} While for a* > 0.24 
one has < < the order is reversed for a* < —0.75. 



III. KINETIC MODEL DESCRIPTION 



The description in the previous Section is based on the Boltzmann equation. It has, however, two shortcomings. 
On the one hand, it is restricted to Maxwell molecules. On the other hand, even for Maxwell molecules, the explicit 
expression for the velocity distribution function is not known. Both limitations are overcome if one resorts to a 
description less detailed than that offered by the Boltzmann equation and employs a model kinetic equation based 
on it. The simplest and best known model kinetic equation is the one proposed by Bhatnagar, Gross, and Krook 
(BGK) p3|. It consists of replacing the true collision operator J[f, /] by a single-time relaxation term of the form 
—v{f — Jl), where v is an effective collision frequency and 



d/2 



mV 



h = n {2^r) e ^{-2k^) (30) 

is the local equilibrium distribution function. The collision frequency v is also a functional of / through its dependence 
on the density and the temperature. While the dependence on n is always linear, its dependence on T varies according 
to the interaction potential under consideration. For instance, in the case of repulsive potentials of the form <p(r) ~ r~^, 
we simply have v cx nT 1 with 7 = | — (d— l)//x p7| . The extreme cases correspond to Maxwell molecules (// = 2(d— 1), 
7 = 0) and hard spheres (fi — > 00, 7 = 4). For this class of repulsive potentials, the Boltzmann equation (Q) is modeled 
as 

i? " a W x Vx ) ^ -MT/nnJ- Jl), Jl = (no/n)/ L , (31) 

where, as before, the subscript denotes initial values. Since the BGK model contains a single parameter (y), it is 
unable to reproduce the correct Boltzmann values of both the shear viscosity an the thermal conductivity coefficients 
simultaneously. In our problem, however, only the shear viscosity is relevant and thus the effective collision frequency 
v can be chosen as v — p/rj-^s, so that the exact NS viscosity is recovered. 



A. Maxwell molecules 



In this subsection we specialize to Maxwell molecules (7 = 0). In that case, the evolution equations for the elements 
of the pressure tensor are identical to those already derived from the Boltzmann equation, Eqs. (||) and (|io|), provided 
that the BGK collision frequency i/q is identified with the one arising from the Boltzmann equation. As a consequence, 
the nonlinear viscosity is again given by Eqs. (|l^) and ([l2]). However, velocity moments of degree higher than two 
no longer coincide in both descriptions. The main advantage of the BGK equation is that it lends itself to an exact 
solution at the level of the distribution function. This solution is expected to provide a fair description of the true 
distribution at least for velocities smaller than or of the order of the thermal velocity (2UbT /m) 1 ' 2 . 



6 



The general solution of Eq. ( |3l| ) for Maxwell molecules is 



7(V,r) = e- VoT e aTdv * v *f(V,0) 

+u f dT'e^ o(T - T V (T - T ' )s ^ y */ L (V,T'), 
Jo 



where the action of the operator exp(ardv x V x ) is 



F{V X ) = e aT F(e aT V x ). 



(32) 



(33) 



We will focus on the long-time distribution function, which becomes independent of the choice of the initial distribution 
/(V, 0) = /(V, 0). To that end, it is convenient to work with the reduced quantities (p0|). The BGK counterpart of 
Eq. @ is 



-a -t^t^x 



A(o*) d 
2 di 



(34) 



whose solution is 



$(€;o*)=7r- 



-d/2 



ds exp 



(35) 



Of course, the same result is obtained by taking the limit vqt — ► oo in Eq. (p2|). Equation (]3q) also allows us to get 



the moments (£ 2fcl ^ 2 ) ■ A simple calculation yields 



r(*i + i)r(fc 2 + 4=i 



r(i)r( d - 



+ — ) 

jt — 2 — - [1 + fci(2a* — A) - fcaA]" 1 



(36) 



if 1 + fci(2a* — A) — feA > 0, being oo otherwise. The sign of 2a* — A = (d — 1)A/(1 — dA) is the same as that of 
a*. Thus the moment (£ 2 ' Ci £_l 2 ) diverges when k\ is sufficiently low and k% is sufficiently high in the case a* > 0, 
while it diverges when k\ is sufficiently high and ki is sufficiently low in the opposite case a* < 0. More specifically, 
(£l fe ) -> oo for a* > if jfe > 1/A(a*) > d; analogously, -> oo for a* < if jfe > l/|2a* - A| > - 1). In 

the particular case of the fourth-degree moments, anc ^ ($1) remain finite but diverges if A is equal to or 

smaller than a critical value A c = — 2), i.e., if a* < a* = —d/A(d— 2). This behavior is reminiscent of the 
one observed in Sec. || from the Boltzmann equation, although there are two main differences: (i) in the case of the 
Boltzmann equation the three fourth-degree moments (i.e., not only diverge and (ii) that happens for a larger 
departure from equilibrium (a* ~ —1.599 versus cC_= —0.75 for d — 3). The dependence of the three moments on a*, 
as predicted by the BGK model, is shown in Fig. ^ It can be observed a good agreement with the results obtained 
from the Boltzmann equation in the region a* > —0.3, especially in the case of For longitudinal rates a* < —0.3, 
however, the deviations become important. 

It is remarkable that in this viscous longitudinal problem the BGK model is able to capture, at least at a qualitative 
level, the existence of diverging moments for values of \a*\ sufficiently large. In the case of the uniform shear flow, 
however, all the moments predicted by the BGK equation are finite pq , pq| , in contrast to the scenario arising from the 
Boltzmann equation p5}-|23]| . The origin of diverging moments can be traced back to the existence of a high- velocity 
tail in the distribution function. To clarify this point, let us consider the two marginal distribution functions 



-1/2 



ds exp 



(37) 



/OO 
-oo 

ds exp 



-(d-l)/2 



(38) 



where £j_ = £ — £ a x is the transverse velocity vector. A simple change of variable in Eq. (p7[) gives 
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ft *\ -1/2^(0,*, £ x ) >-2/3|| a , *s _ 1 1 , on , 

^"^ ;fl)=7r ' ^l (a)= 2-2^A' (39) 



where 

r 

Here, 



Ma* r (Ah^), «*>0, f40 s 

|| (« jr(/3||)-r(/?||,£2), a*<0. (4 ° j 



r(/3,x) = / dyy^e^ (41) 



is the incomplete gamma function p?]. Analogously, 



ipx ^ x; a*)=n-^ ^^ ^, /3 ± (a*) = ^ + i (42) 



where 



p r„* , 2 ^_ fr(^)-r(/? ± ,ei), a* >o, 



(43) 



From Eqs. (||) and ( |4p| ) it follows that <pn ~ £ x 2/3,1 in the limit £| — * oo if a* < 0, and so (£| fc ) — >ooiffc>/3ii — 4 = 



2 

|2a* - A| _1 . Similarly, tp ± ~ f" 2 ^ in the limit -> oo if a* > and then _ MX jiffc> J g L _<^l = A" 1 . 

It is interesting to note that the exponents /3m and /3_l remain finite in the limit of infinite \a*\: liniQ.^oo (3±(a*) = 
(<f-l)lim«._^_ ooJ 8||(a*) = (3d-l)/2 and lima-^-oo 0x(a*) = (d-^lim,,.^ ^(a*) = (<f-l)/2. 

Apart from an algebraic high- velocity tail, the distribution function may exhibit a singular behavior in the opposite 
limit of vanishing velocities. Equation (|35| ) shows that lim^o "£(£; a *) = oo if 1 — a* + |A(a*) < 0, which corresponds 
to a* < — [y/d(3d — 2) — d+2]/2(d— 1) and a* > \y/d(3d — 2) + d— 2]/2(d— 1). A similar phenomenon of overpopulation 
of "rest" particles occurs in the uniform shear flow state P5fl . Again, it is useful to consider the marginal distributions 
in the analysis of this effect. In the case of <p\\ the divergence happens when both /3y and a* are positive, i.e., for 
a* > 3d/(3d — 1), while in the case of ip± the singular behavior takes place when (3± > and a* < 0, i.e., for 
a* < -d(d+ l)/{d - l)(3d- 1). By using the properties @|t| 




lim T(P,x) = { -p-^, f3<0, (44) 

one gets from Eqs. (^) and ( fl2"| ) the following asymptotic behaviors: 

r ^(i^^ + iA)- 1 , a * 

Jm^;a*) = \ «-)'im;G 2P \ a* > (45) 



a 



3d 
3d-l ' 



lira = vr-^ 1 )/ 2 ^^^, a* < 'J^ry (46) 

Figure|| shows the ratios with respect to local equilibrium i?y(^ x ;a*) = (pi\ (£ x ; a*)/ip\\ (£ x ; 0) and R±(£ ± ;a*) = 
95_l(^j_; a*)/i/?j_(^j_; 0) for a* = —1 and a* = 1.5 in the three-dimensional case. It can be observed that at a* = — 1 
(a* = 1.5) the function cpn (<p±) develops a high-velocity tail, while the function tp± (cpn) diverges as the velocity 
vanishes. 
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B. Repulsive potentials. Hard spheres 



Now we consider more general repulsive potentials characterized by 7 > 0. The BGK model for our problem is 
given by Eq. (|3l|). Its general solution is 



/(V,r) = e - s Me a *- v '/(V,0) 



-vq 



f At' 


\T{r')] 


la 


. To . 



where 



a _[ s ( T )_ s ( T ')] a(T-r')dv K V x 



/z(vy), 



(47) 



s(r) = i/ 



dr' 


\T(t')] 


'0 





(48) 



is the number of collisions per particle. In the limit 7^0, s(r) = v^t and Eq. ( f47j ) reduces to Eq. (|3^). On the other 
hand, for 7 > Eq. (^) is not closed since it requires the knowledge of the time dependence of the temperature. 
By multiplying both sides of Eq. ( [47| ) by V 2 and integrating over velocity one can get a closed integral equation for 
T(r). Such an equation is however quite involved and then it is more transparent to work directly with the evolution 
equation itself, Eq. ( pl| ) . From this equation it is straightforward to find that the evolution of p and P xx is still given 
by Eqs. (^) and (flO|), except that now uq is replaced by a time-dependent collision frequency vo[p(t) / po] 7 • Therefore, 
the pressure p — n^ksT obeys a nonlinear second-order differential equation 



dr 2 



2a 



vq\ — 



d 

—p + -av I — 
or d \po 



P, 



(49) 



subject to the initial conditions p(0) = po, dp/dr\ T=0 = —(2a/d)P xx (0). For asymptotically long times (r — > 00), the 
solution of Eq. ( |is| ) behaves as p(t) ~ exp(— 2ar/d) if a < and as p(r) ~ (1 + jvor/d)^ 1 ^ if a > 0. In both cases 
the accumulated number of collisions s(t), Eq. (Eq), goes to infinity as r -> 00. In particular, if a < 0, the typical 
number of collisions per particle during the finite interval < t < |a| _1 becomes infinite, a property already seen in 
Sec. |0] for Maxwell molecules. 

In order to get the nonlinear viscosity as a function of the longitudinal deformation rate, we must work with the 
reduced rate a* = Ajv = a/^o(p/po) 7 rather than with the time variables r or t. Note that lim r ^ 00 a*(r) = 00 if 
a > 0, while lim,-^^ a* (r) = if a < 0. With this change of variable one has 



da 



~P 



P 

7a* 



(50) 



—P xx = _*P_ 

da* xx -fa* 



1 



1 

2a* 



1 - 



P 
p 

1 XX 



(51) 



This gives rise to the following ordinary differential equation for the reduced nonlinear viscosity defined in Eq. (|l4|): 

.d 



2 ia * 2 [ 1 - 2^— -a*rf 



di]* 
da* 



4— r^(l - 7)a*V 2 + [d + 2(d - 2 + 7)0*] rf - d = 0. 



(52) 



For small a* the solution is rf w 1 — [2(d — 2 + j)/d]a*, while the asymptotic behavior of -q* for large |o*| is 
n* fa [d/2{d-l)]a*- 1 (l-a*- 1 /2) for a* > and n* w (d/2)|a*|^ 1 [l - |a*|- 1 /2(l + 7)] for a* < 0. The leading terms 
are independent of 7 and correspond to lmi a * ^ +OQ P xx /p = and \im a * P xx /p — d, respectively. Interestingly 
enough, Eq. (p2|) , particularized to = 3, is equivalent to the one derived by Karlin et al. || from their invariance 
principle and Grad's method. In other words, the invariance principle under the microscopic and macroscopic dynamics 
is an approximation that, at least in this problem, yields the same nonlinear viscosity as the one predicted by the 
BGK model. The latter approach, nevertheless, has the advantages of being conceptually simpler and providing the 
full velocity distribution function. 

It must be noted that Eq. (|52|) possesses as many solutions as particular initial conditions. Each particular solution 
is specified by assigning a given value of the viscosity rf at the initial (reduced) longitudinal rate a^ = ajvQ. This 
situation is analogous to the one discussed in Ref. |2q1 for the uniform shear flow case. Since the irreversible time 
evolution of the system leads to a monotonic increase of a*, Eq. (p2]) must be solved for a* > at. For positive rates, 
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this implies the range < Oq < a* , but for negative rates the range is Oq < a* < 0. All the particular solutions, 
however, tend towards a special solution (the hydrodynamic one) for sufficiently long times, i.e., for a* ^> a,g if Uq > 
and for \a*\ <§C |ag| if Oq < 0. This special solution representing the hydrodynamic or normal reg ime can be identified 
in principle by the Chapman-Enskog series (|^|). Insertion into the differential equation ( |52|) yields the following 
recurrence formula, 

d — 2 + 717 d — 1 \ 

Cn = -2 Cn-1 - 4 rf2 2 , CmCn-2-m [1 - (n - 1 - TOj 7 J . (53) 

' m=0 

This equation shows that the coefficient c„ is a polynomial in 7 of degree n. The first four coefficients in the case of 
a three-dimensional system of hard spheres (d = 3, 7 = 5) are ci = — 1, c 2 = §,03 = — 1|, and C4 = Further 
computation of the coefficients shows that, except in the case of Maxwell molecules (7 = 0), the expansion jl^ ) 
is only asymptotic. For large n, the magnitude of the coefficients c n grow so rapidly that the second term on the 
right-hand-side of Eq. ( p3| ) can be neglected, so c„/c„_i w —2nj/d. The ratio — c n /c n -i, n = 1-20, is plotted in Fig. 

for three-dimensional hard spheres. The linear growth of the ratio is already apparent for n > 4. The divergence of 
the Chapman-Enskog expansion for 7 > also takes place in the uniform shear flow problem ]25| , |2"8| ] . 

Since the series (|l£]) is only useful if truncated and applied to small a* , a different strategy is needed to get the 
hydrodynamic rf for finite a*. One possibility is to expand 77* around the point at infinity (namely in powers of 
a* -1 ). Such an expansion proved to be convergent in the case of uniform shear flow for shear rates larger than a 
certain value [^5|,^8j. From a practical point of view, however, this method is not very convenient because many terms 
would need to be retained in order to get reliable results in the range of interest ( say \a*\ ~ 1), even if the expansion 
converges. A second possibility is to solve numerically the differential equation (|52j ) with the boundary conditions 
lim a *^ + 77*(oq) = 1 (for a* > 0) and lim a »_>_oo 77* (ap) = (f° r a * < 0)- On the other hand, it seems more convenient 
to follow the approach proposed in Ref. Jl!|, which consists of representing the nonlinear viscosity as an expansion in 
powers of the interaction parameter 7: 

N 

n*(a*) = km V ^(a*), V (N \a ) = V r,„(a*) 7 ", (54) 



I 



where ?7o( a *) is the nonlinear viscosity for Maxwell molecules, Eq. (J15|). Since the Chapman-Enskog coefficients c„ 
are polynomials in 7, the series (|54|) can be interpreted as a rearrangement of the series (16). In other words, if we 
write 



m=0 

then 

oo 

Vn(a*) = ^2 c mn a* m . (56) 

m—n 

As a consequence, the truncated series rf N \a*) is exact through order a* N , i.e., 77* (a*) — ri^ N \a*) = 0(a* N+1 ). 
Insertion of Eq. (M) into Eq. ([32]) yields the following recurrence formula 



i]n(a*) 



2a 



d+2(d-2)a* +8^a* 2 r l0 (a*) 
n-l 



1 - 2^— -a* m (a*) 



[r) n -i(a*) + a*ri'n-i(a*)] 



~ 2 ~~[~ a * X! ^mia*) [vn-i- m (a*) + a*J7^_ x _ m (a*) - r?„_ m (a*)] \ , (57) 

m— 1 ) 

where the prime denotes a derivative with respect to a*. The coefficient r/i for d = 3 was the only one considered 
by Karlin et al. P, |Tl| , p9| . The coefficients r) n for n = 0, 1, 3, 5, 6 and d = 3 are plotted in Fig. Up to n — 3 the 
coefficients remain small, but the magnitude of 775 and, especially, that of rje reach rather high values, thus suggesting 
the asymptotic character of the expansion (p4), at least for negative a*. This is not surprising if one takes into 
account that, while all the truncated series rf^^(a*) are regular at a* — 0, the full viscosity 17* (a*) is singular at 
a* = 0. Notwithstanding this, since the maximum value of 7 is 7 = \, it turns out that the functions rj^ N \a*) 
with iV = 3 or N = 4 can be considered as rather good approximations of 77* (a*). This is quite apparent in Fig. @, 
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where j/ 3 )(a*) and ri^(a*) practically overlap in the case 7=3 (corresponding to a repulsive potential with fi = 12) 
and are hardly distinguishable in the case of 7 = \ (hard spheres) . Figure [5] also shows that the nonlinear viscosity 
is almost insensitive to the interaction model in the case of an expansion (a* > 0). On the other hand, when the 
physical situation corresponds to a condensation of the gas (a* < 0), the hardness of the repulsion plays a relevant 
role, especially around the maximum (a* « —0.4) p0| . 

Before closing this Section, it is worthwhile noting that Uribe and Garci'a-Cohn [|l2| used Grad's (nonlinear) moment 
method to get an expression for 77* (a*) that is dramatically at odds with the results obtained in this paper for a* < 0. 
According to their results, 77* (a*) monotonically increases as a* becomes more and more negative and finally reaches 
a plateau 77* — * 49 in the limit a* — > — 00. However, these results are strongly inconsistent with the physical condition 
P = [Pxx + (d- l)P yy \/d > P xx /d, which implies [cf. Eq. (|T|)] that 77* (a*) < (d/2)\a*\- 1 if a* < 0. Since rj* has an 
upper bound that goes to zero in the limit a* — > —00, then limQt^-oo r/*(a*) = necessarily. On the other hand, from 
Eqs. (45) or (47]of Ref. |H it follows that p < P xx /3 ( i.e, Pyy < 0) if a* < — These inconsistencies of the results 
derived in Ref. |12| ] are likely associated with the assumption of a stationary situation in this compressible flow. 

IV. CONCLUSIONS 

This paper has dealt with a simple viscous longitudinal flow characterized by an unsteady velocity profile u x (x, t) = 
ax/(l + at) and a uniform density n(t) = no/(l + at). The situation with a > corresponds to an expansion of the gas, 
while the case a < describes a condensation phenomenon. By using kinetic theory tools (Boltzmann equation and 
BGK kinetic model), exact results have been derived for the generalized or nonlinear viscosity, the velocity moments, 
and the velocity distribution function. The following points summarize the main conclusions of the present study. 

• By an adequate change of velocity and time variables, the problem is seen to be formally equivalent to that of 
a uniform gas with a steady density, in which the particles are under the action of a longitudinal driving force 
F = — maVxX.. According to this viewpoint, the particles are decelerated or accelerated along the longitudinal 
direction, depending on the sign of a. As a consequence, the temperature monotonically decreases in time if 
a > 0, while it increases if a < 0. 

• The relative difference between the normal stress P xx and the hydrostatic pressure p for long collision times is 
characterized by a (dimcnsionless) viscosity coefficient rf (a* ) , which is a nonlinear function of the longitudinal 
deformation rate relative to an effective collision frequency. The (Chapman-Enskog) expansion of 77* in powers of 
a* is, in general, only asymptotic. An exception is provided by Maxwell molecules, in which case the expansion 
converges for \a*\ < \. 

• A thinning effect is present for a* > 0, i.e., 77* monotonically decreases as a* increases. For a* < 0, however, 
77* starts increasing with \a*\ (thickening effect), reaches a maximum (at a* = — i for Maxwell molecules and 
around a* ~ —0.4 for hard spheres), and then decreases for more negative longitudinal rates. 

• In the case of an expansion (a* > 0), the nonlinear viscosity 77* is practically "universal". On the other hand, 
its behavior in the case of condensation (a* < 0) is rather sensitive to the interaction potential. In particular, 
the harder the potential the higher the maximum value of 77* (for instance, 77j^ ax = 1.125 for Maxwell molecules 
and 77^ ax ~ 1.46 for hard spheres). 

• The results for 77* derived from the Boltzmann equation for Maxwell molecules and from the BGK model for 
more general potentials coincide with those derived by Karlin el al. || from Grad's method and the application 
of their invariance principle under microscopic and macroscopic dynamics. It would be interesting to explore 
whether such an equivalence extends to the similar but more complicated problem of uniform shear flow as well. 

• The shape of the velocity distribution function (for Maxwell molecules) has also been analyzed by scaling 
the velocities with the (unsteady) thermal velocity. This gives rise to a new term in the kinetic equation that 
represents a non-conservative thermostat force that cancels the cooling (a* > 0) or heating (a* < 0) produced by 
the driving force. This exact equivalence between the free system and the thermostatted one is analogous to that 
taking place in the uniform shear flow and is restricted to Maxwell molecules (collision frequency independent 
of the velocity). For other interactions the equivalence is only approximate, but yet the thermostatted problem 
is worth studying by itself. It would be quite interesting to carry out nonequilibrium molecular dynamics 
simulations of hard spheres subject to the simultaneous action of the driving and thermostat forces, in order to 
measure the nonlinear viscosity and related phenomena in dense gases. This would complement the extensive 
simulation studies of the uniform shear flow. 
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• The exact fourth-degree (scaled) velocity moments derived from the Boltzmann equation diverge for sufficiently 
negative values of the longitudinal rate (a* < —1.6). This indicates the existence of an algebraic high- velocity 
tail, especially for the longitudinal component of the velocity, for negative a* . 

• The above singular behavior of the moments is described, at least qualitatively, by the exact solution of the BGK 
model. This is a very remarkable feature, since a similar behavior showing up in the uniform shear flow was not 
captured by the kinetic model. Also, note that Grad's method is unable to predict an algebraic high-velocity 
tail, as it approximates the distribution function by a Gaussian times a polynomial. The BGK solution obtained 
in this paper predicts an algebraic tail in the marginal distribution of longitudinal velocities if a* < and a 
weaker tail in the marginal distribution of transverse velocities if a* > 0. The latter tail implies that, while all 
the moments of degrees equal to or smaller than six (for a three-dimensional system) are finite when a* > 0, the 
moments of eighth degree diverge if a* > 1.125. The investigation of whether or not this prediction is confirmed 
by the Boltzmann equation will be the subject of a separate paper. 

• The explicit expression of the BGK velocity distribution function allows one to unveil a different type of singular 
behavior that does not have, however, a direct influence on the velocity moments since it is associated with the 
limit of small velocities. More specifically, the distribution of vanishing velocities diverges if the longitudinal rate 
is sufficiently positive (a* > 1.40) or sufficiently negative (a* < —0.90). This effect is also present in the marginal 
distribution of longitudinal velocities (for a* > 1.125) and in that of transverse velocities (for a* < —0.75). 
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a 

FIG. 1. Plot of the reduced fourth-degree moments (k — 0), (£^£x) = 1)) ar >d (£x) (k = 2) as functions of the reduced 
longitudinal rate a* in a three-dimensional system of Maxwell molecules. Solid lines refer to the exact results derived from the 
Boltzmann equation, while dashed lines are predictions of the BGK kinetic model. 
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£ 2 £ 2 

FIG. 2. Marginal velocity distributions normalized with respect to local equilibrium, R\\ (solid lines) and R±{£±) (dashed 
lines), for a* = —1 and a* = 1.5, as predicted by the BGK kinetic model for a three-dimensional system of Maxwell molecules. 
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n 

FIG. 3. Ratio -c n /c n _i between successive coefficients in the Chapman-Enskog expansion of the nonlinear viscosity for a 
three-dimensional system of hard spheres, according to the BGK kinetic model. 
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a 

FIG. 5. Nonlinear viscosity for three-dimensional systems of Maxwell molecules (7 = 0), particles interacting via an 
r~ 12 -potential (7 = |), and hard spheres (7 = |). The solid line is the exact result derived from the Boltzmann equa- 
tion and the BGK model, while the dashed and dotted lines are the approximations rf z > and rf 4 ', respectively, as obtained 
from the BGK kinetic model. 
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